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Abstract 

Gravity stands out among the fundamental interactions because of 
its apparent incompatibility with having a quantum description. More- 
over, thermodynamic aspects of gravitation theory appears as puzzling 
features of some classical solutions such as black holes. These and other 
aspects of gravitational theories have recently lead to the proposal that 
gravity might not be a fundamental interaction but rather an emergent 
phenomenon, a sort of hydrodynamic limit of some more fundamental the- 
ory. In order to further explore this possibility we shall here discuss two 
systems where such emergence of a gravitational dynamics is observed. 
We shall consider first the case of a non-relativistic Bose- Einstein conden- 
sate and then a more abstract model based on scalar fields living on a 
Riemannian manifold. This will allow us to put in evidence the general 
issues related to emergent gravity scenarios with a particular attention to 
the role and nature of Lorentz and diffeomorphism invariance. 



1 Introduction 

In spite of being the first force of Nature to be understood in physical terms, 
gravity is somehow still a riddle for physicists. Not only it keeps evading a full 
quantum description as well as any form of unification with the other interac- 
tions, it also puzzles us with profound questions and unexpected features. We 
will not attempt here to present a complete list of these startling aspects of 
gravitation theory, but we can recall, for example, the surprising connection be- 
tween gravity and thermodynamics associated to black hole physics [1] |^ [3] as 
well as the deep questions associated to the nature of inertia and time [H O [S] . 

In recent years, a new approach to these old problems has been gaining mo- 
mentum and many authors have been advancing the idea that gravity could 
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all in all be an intrinsically classic/large scale phenomenon similar to a con- 
densed matter state made of many atoms [7]. In this sense gravity would not 
be a fundamental interaction but rather a large scale/number effect, something 
emergent from a quite different dynamics of some elementary quantum objects. 
In this sense, many examples can be brought up, starting from the causal set 
proposal 01 passing to group field theory [3] or the recent quantum graphity 
models [10] and other approaches (see e.g. [TT]). 

All these models and many others share a common scheme: they consider a 
fundamental theory which is not General Relativity and examine, using differ- 
ent techniques often borrowed from condensed matter physics, how space, time 
and their dynamics could emerge in some regime. It is perhaps important to 
remark that such "emergent gravity" scenarios should not be seen as alterna- 
tives to quantum gravity proposals (e.g. superstrings theory or loop quantum 
gravity) rather one should think of them as different incarnations of a more 
general paradigm about how classical spacetime geometry and dynamics could 
be recovered from such quantum gravity scenarios. 

In this sense a leading inspirational role also been played by a parallel stream 
of research which goes under the name of "analogue models of gravity" [12]. 
These are condensed matter systems which have provided toy models showing 
how at least the concept of a pseudo-Riemannian metric and Lorentz invariance 
of matter equations of motion can be emergent. For example, non-relativistic 
systems which admit some hydrodynamics description can be shown to have 
perturbations (phonons) whose propagation is described, at low energies, by 
hyperbolic wave equations on an effective Lorentzian geometry [12j . While 
these models have not provided so far also an analogy of emergent gravitational 
dynamics equations they do have provided a new stream of ideas about many 
other pressing problems in gravitation theory (see for example recent works on 
the origin of the cosmological constant in emergent gravity |13|). The stud- 
ies presented here should then be considered as exploratory toy models aimed 
at gaining an understanding of the crucial ingredients necessary for a general 
emergence paradigm to work. 

2 Analogue models as a test-field for emergent 
gravity 

Analogue models for gravity have provided a powerful tool for testing (at least in 
principle) kinematical features of classical and quantum field theories in curved 
spacetimes [12]. The typical setting is the one of sound waves propagating 
in a perfect fluid [TU [15]. Under certain conditions, their equation can be 
put in the form of a Klein-Gordon equation for a massless particle in curved 
spacetime, whose geometry is specified by the acoustic metric. Among the 
various condensed matter systems so far considered, Bose-Einstein condensate 
(BEG) [ISl dZ] had in recent years a prominent role for their simplicity as 
well as for the high degree of sophistication achieved by current experiments. 
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In a BEC system one can consider explicitly the quantum field theory of the 
quasi-particles (or phonons), the massless excitations over the condensate state, 
propagating over the condensate as the analogue of a quantum field theory of 
a scalar field propagating over a curved effective spacetime described by the 
acoustic metric. It provides therefore a natural framework to explore different 
aspects of quantum field theory in various interesting curved backgrounds (for 
example quantum aspects of black hole physics [TH [19] or the analogue of the 
creation of cosmological perturbations [201 HI])- Unfortunately, up to now, the 
analogy with gravity is only partial: there is no analogy with some sort of 
(semiclassical) Einstein equations, since it has not been possible to put the fiuid 
equations, which are those describing the dynamics of the acoustic metric, in a 
geometrical form which could eventually lead to a complete dynamical analogy 
with general relativity [22]. Our first task here is to show how to fill this gap 
in the case of BEC and to gather, from the understanding of the emergence of 
a BEC analogue gravitational dynamics, general lessons about possible features 
of "emergent gravity" scenarios. 

3 Emergent spacetime in BEC: a review 

In BEC, the effective emerging metric depends on the properties of the con- 
densate wave- function. One can expect therefore the gravitational degrees of 
freedom to be encoded in the variables describing the condensate wave-function 
[17], which is solution of the well known Bogoliubov-de Gennes (BdG) equation. 
The dynamics of gravitational degrees of freedom should then be inferred from 
this equation, which is essentially non-relativistic. The gravitodynamics of the 
BEC should therefore be the analogue of some sort of Newtonian gravity, and 
we shall reinterpret the BdG equation as a modified Poisson equation. 

The "emerging matter" , the quasi-particles, in the standard BEC, are phonons, 
i.e. massless excitations. Since we expect the quasi-particles to be the matter 
source in the Poisson equation, we run a priori into a problem: massless particles 
are not treatable in the framework of Newtonian mechanics. To avoid this issue, 
we shall then introduce a new term in the BEC Hamiltonian which will softly 
break the usual U{1) symmetry and therefore will allow the quasi-particles to 
acquire mass. 

3.1 The hydrodynamic limit 

Let us start by very briefly reviewing the derivation of the acoustic metric for 
a BEC system, and show that the equations for the phonons of the condensate 
closely mimic the dynamics of a scalar field in a curved spacetime. In the dilute 
gas approximation, one can describe a Bose gas through a quantum field ^' 
satisfying 
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m is the mass of the atoms, a is the scattering length for the atoms and k 
parameterises the strength of the interactions between the different bosons in 
the gas. It can be re-expressed in terms of the scattering length a as 

K(a) = (2) 

As usual, the quantum field can be separated into a macroscopic (classical) 
condensate and a fluctuation: = ijj + (p, with {^) = Then, by adopting 
the self-consistent mean field approximation 

0^(p(p ~ 2{0^(p) (p + {(pip) (p\ (3) 

one can arrive at the set of coupled equations: 

ifi— V(t,x) = ( -— + Vcxt(x) -h/cnJ V(i,x) 

-hK{2nV(t,x)-hmV'*(i,x)}; (4) 

. f h^ 

+K rriT (p\t,-x.). (5) 

Here 

ric = |V'(i,x)i^; mc = V^(i,x); (6) 

n={(p^V)\ fh={!f(p); (7) 

riT = nc + n; tut = mc + rh. (8) 

In general one will have to solve both equations for ip and (p simultaneously. 

The equation for the condensate wave function ijj is closed only when the back- 
reaction effects due to the fluctuations are neglected. (The back-reaction being 
hidden in the quantities rh and h.) This approximation leads then to the so- 
called Gross Pitaevskii equation. 

Adopting the Madelung representation for the wave function ijj of the con- 
densate 

tp{t,x) = ^/nc{t,x) ex.p[-ie{t,x)/h], (9) 

and defining an irrotational "velocity field" by v = 'V9/m,, the Gross- Pitaevskii 
equation can be rewritten as a continuity equation plus an Euler equation: 

d 

— Tic + V • (nev) = 0, (10) 
+ yext(i,x)+«nc-— ^^1=0. (11) 
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These equations are completely equivalent to those of an irrotational and inviscid 
fluid apart from the existence of the so-called quantum potential 



^quantum = -ft V ^71^/(2™^/^^) 

which has the dimensions of an energy. Note that 



n. V,; 



2m 



4m 



which justifies the introduction of the so-called quantum stress tensor 

re 



quantum 



Am 



(12) 



(13) 



(14) 



This tensor has the dimensions of pressure, and may be viewed as an intrinsically 
quantum anisotropic pressure contributing to the Euler equation. If we write 
the mass density of the Madelung fluid as p = m ric, and use the fact that the 
flow is irrotational then the Euler equation takes the form 



P 



V -f (v • V)v 



-l-p V 




+ V 


r 9 1 




m 




2m2 



y . ^quantum ^ (15) 



Note that the term V^xt/^i has the dimensions of specific enthalpy, while k/9^/(2to) 
represents a bulk pressure. When the gradients in the density of the condensate 
are small one can neglect the quantum stress term leading to the standard hy- 
drodynamic approximation. Because the flow is irrotational, the Euler equation 
is often more conveniently written in Hamilton-Jacobi form: 



4' 



2m 



-I- yoxt(i,x) -I- uric 



2m 



= 0. 



(16) 



Apart from the wave function of the condensate itself, we also have to account 
for the (typically small) quantum perturbations of the system ([5]) . 

Let us consider now the quantum perturbations above the condensate. These 
can be described in several different ways, here we are interested in the "quantum 
acoustic representation" 



(p{t, x) 



-w/h 



1 



ni 



(17) 



where ni, 0i are real quantum fields. By using this representation Equation ^ 
can be rewritten as 



dtui 



1 



V ■ [niVe + ncVei] = 0, 



1 ft2 

dtOi H Ve ■ V9i + K(a) ni - — D2ni 

m 2m 



(18) 
(19) 
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Here D2 represents a second-order differential operator obtained from linearizing 
the quantum potential. Explicitly: 



2 " 



,+1/2 



)] m 



-1/2 



n-i) 



(20) 



The equations we have just written can be obtained easily by linearizing the 
Gross-Pitaevskii equation around a classical solution: ric — > ric + rii, </) ^ 
+ 01. It is important to realise that in those equations the back-reaction of the 
quantum fluctuations on the background solution has been assumed negligible. 
We also see in Equations (fT8l[T9|) . that time variations of Vcxt and time variations 
of the scattering length a appear to act in very different ways. Whereas the 
external potential only influences the background Equation (fT6|) (and hence 
the acoustic metric in the analogue description), the scattering length directly 
influences both the perturbation and background equations. From the previous 
equations for the linearised perturbations it is possible to derive a wave equation 
for 01 (or alternatively, for rii). All we need is to substitute in Equation (|18p 
the ni obtained from Equation p9)) . This leads to a PDE that is second-order 
in time derivatives but infinite order in space derivatives - to simplify things we 
can construct the symmetric 4x4 matrix 



f 



iO 



(21) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Then, intro- 
ducing (3-|-l)-dimensional space-time coordinates 



= (t; x') 

the wave equation for di is easily rewritten as 

df^ir-" dji) = 0. 

Where the f^'^ are differential operators acting on space only: 



TO 

-1 



r = 




K(a) - 


ft2 

2m 






K{a) - 


2m 


r = 


m 


K{a) 




m 


TO 



2m 



K{a) 



2m 



Do 



(22) 
(23) 

(24) 
(25) 
(26) 
(27) 



Now, if we make a spectral decomposition of the field 9i we can see that for 
wavelengths larger than ^ = h/mcsound (C corresponds to the "healing length" , 
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as wc will explain below and Csound(a, nc)^ = -^^^7)^^), the terms coming from the 
linearization of the quantum potential (the D2) can be neglected in the previous 
expressions, in which case the can be approximated by scalars, instead of 
differential operators. (This is the heart of the acoustic approximation.) Then, 
by identifying 



-9 9" 



(28) 



the equation for the field 61 becomes that of a (massless minimally coupled) 
quantum scalar field over a curved background 



A0i 



1 



-g g'^" d.) 01 = 0, 



(29) 



with an effective metric of the form 



9ixu{t,-ii) = 



m Csound(a,nc) 



-{Csound(a,'^c)^ - v"^} 



-Vi 



Si 



(30) 



Here the magnitude Csound{nc, a) represents the speed of the phonons in the 
medium: 

Csound (a, ric) = — — — , (31) 
and Vi is the velocity field of the fluid flow, 

1 



m 



Vi 



m 



(32) 



3.2 Lorentz violation in BEC 

It is interesting to consider the case in which the above "hydrodynamical" ap- 
proximation for BECs does not hold. In order to explore a regime where the con- 
tribution of the quantum potential cannot be neglected we can use the so called 
eikonal approximation, a high-momentum approximation where the phase fluc- 
tuation 9i is itself treated as a slowly- varying amplitude times a rapidly varying 
phase. This phase will be taken to be the same for both ni and 61 fluctuations. 
In fact, if one discards the unphysical possibility that the respective phases dif- 
fer by a time varying quantity, any time-independent difference can be safely 
reabsorbed in the definition of the (complex) amplitudes Ao, Ap. Specifically, 
we shall write 

?i(t,x) = Re{Ae exp{-i(j))} , (33) 
ni(i,x) = Re {Ap exp{-i(j))} . (34) 

As a consequence of our starting assumptions, gradients of the amplitude, and 
gradients of the background fields, are systematically ignored relative to gradi- 
ents of (/). (Warning: What we are doing here is not quite a "standard" eikonal 
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approximation, in the sense that it is not appHed directly on the fluctuations of 
the field ipitjX) but separately on their amplitudes and phases pi and We 
adopt the notation 



Then the operator D2 can be approximated as 
L»2 ni = 



(35) 

(36) 
(37) 
(38) 



A similar result holds for D2 acting on Q\. That is, under the eikonal approxi- 
mation we effectively replace the operator D2 by the function 



1 



D2 ^ -l^n-'k 

For the matrix /^^ this effectively results in the replacement 



^00 



K{a) 
K(a) 



4m ric 



4m n, 
K{a) + 



4m ric 



4m ric 



(39) 

(40) 
(41) 
(42) 
(43) 



(As desired, this has the net effect of making f^'^ a matrix of numbers, not 
operators.) The physical wave equation (|23p now becomes a nonlinear dispersion 
relation 

+ if' + /«) CO h + P h k, = 0. (44) 

After substituting the approximate D2 into this dispersion relation and rear- 
ranging, we see (remember: — ||fc|p = ki kj) 



That is (with vi = — ViOo) 



Urk^ 



K(a) 



4mnc 



K hf = 0. 



(w - vl kiY 



Urk^ 



n{a) + 



4mnc 



(45) 



(46) 
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Introducing the speed of sound Cgound this takes the form: 

Having described the effective metric in BEC and its hmit of appUcabihty 
we are now in hand, the analogy is fuUy estabhshed, and one is now in a position 
to start asking more specific physics questions. 

4 Emergent gravity in BEC 

The program of extracting some sort of Poisson equation out of the Bogoliubov- 
de Gennes formahsm cannot be carried on in a standard BEC. Indeed, phonons 
are massless excitations, and hence, since we want them to enter the Poisson 
equation as a source term, we have to circumvent the impossibihty of treating 
massless particles in the framework of Newtonian mechanics. This is easily done 
by making phonons massive. Concretely, this is done by introducing a new term 
in the Hamiltonian which will softly break the usual U{1) symmetry associated 
to number conservation and therefore will allow the quasi-particles to acquire a 
mass. Essentially, the quasi-particles will be pseudo-Goldstone bosons [23l [24] : 
their spectrum, instead of being gapless, is gapped. One can expect, then, that 
instead of massless quasi-particles, the collective modes above the condensate 
will be massive. 

In order to do so, the standard Hamiltonian Hq described previously needs 
to be slightly modified, by introducing a term which is (softly) breaking the 
U{1) symmetry in ([T]). 

Ho^H = Ho + Hx, ^^ = -\j '^'^ (^'(x)*(x) -f ^'t(x)*t(x)) . (48) 

The parameter A has the same dimension as /x. With this new Hamiltonian, the 
non-linear equation ([1]) becomes 

d ~ ~ ^ ^ - 

ih—^ = [iJ, m = V^* - + kI^'P* - A*^ (49) 

at 2m 

The addition of this term implies that the whole dynamics of the system 
must be reconsidered. The analysis presented in the subsection (j4.2p will show 
how H\ generates a mass for the quasi-particle. Even though H\ both creates 
and destroys pairs of atoms, it is not difficult to check that H\ is not commuting 
with the number operator N, 

[Hx,N] = j rf^x(*(x)*(x) - *^(x)\^t(x)^ (50) 

while unitarity is preserved. In fact, when applied on a state with a definite 
number of atoms n we have: 

|n) |n-2) + |n + 2), (51) 
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which means that an eigenstate of the number operator evolves into a super- 
position of states with different occupation numbers. However, the expectation 
value of the number of operator on its eigenstates is still constant 

ih—{n\N\n) = {n\[N,H]\n) = {n\[N ,Hx]\n) {n\n - 2) - {n\n + 2) ^0. (52) 

Given this crucial difference with the standard description of BECs, a for- 
malism like the particle-number-conserving one |251 126] cannot be immediately 
used for these new models. However, given the important improvements in the 
description of inhomogeneous condensates provided by this formalism it will be 
interesting to extend the particle-number-conserving method in a suitable way, 
in order to be able to control to which accuracy we can trust the standard mean 
field approximation we are using. In fact, we are assuming that the addition 
of the new term in the Hamiltonian will not destroy the stability properties of 
the mean field theory, and in particular the fact that the mean field theory, 
for nearly homogeneous condensates and trapping potentials, in weakly time- 
dependent regimes, does offer a good description of the condensate dynamics. 

This point of view (and its implications for the condensation mechanism) 
deserves some further specifications. While the condensation process can be 
easily understood as a macroscopic occupation number of an energy level, there 
are several approaches to describe it mathematically. The mean field approach 
is particularly convenient: we say that the system of N bosons has condensed 
whenever the field ^' develops a non-zero vacuum expectation value (vev) 

{n\i!p) = i^, (53) 

where ■0 is the condensate wave-function. If this mean field is non-vanishing, 
we have that the two point correlation function 

G(x,y) = (*^(x)*(y)) « ^*(x)^(y), (54) 

tends to a non-zero constant when x, y are infinitely separated, i.e. the system 
develops long range correlations [57j. 

The mean field method is based on the assumption that the ground state 
of the system is not the vacuum state of the Fock space, |0), but rather it is 
similar to a coherent state. For a single mode Fock space, a coherent state is 
defined to be: 

\z) =e~\'\''^e''^'\{)), (55) 
and it is easy to see that it is an eigenstate of the annihilation operator: 

b\z)=z\z). (56) 

In the case of BEG, the fact that the state uo(x) is macroscopically occupied 
(i.e. there are Nq bosons in the state 0, with Nq/N « 1) can be formalized by 
taking: 

«e-^«/2e-''(^«)''''^J|0). (57) 
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On this states, the field operator ^ behaves hke a c-number: 

« ^ouo{x)\n), (58) 

where the approximation is due to the fact that interactions are introducing 
some corrections. This property of the ground state motivates the sphtting 
of the field operators into the classical part, which deals with the condensed 
phase, and small residual fluctuations, describing the states which are close to 
the ground state: 

^ ^^pi + X. (59) 

This is the basic idea of the mean field approximation: the field operator 5* is 
approximated by a classical field, which is describing the condensate, while the 
fiuctuations, assumed to be small, are still encoded in a field operator. 

The fact that the solutions to the Gross-Pitaevski equation leads to a non- 
vanishing mean field, in light of the discussion about the realization of a regime 
of long range correlations (see ((54|) ). ensures that a condensation has taken 
place. In this sense, the addition of the new term into the Hamiltonian should 
not forbid the condensation, provided that A, /i, k are such that the condensate 
wavefunction can be different from zero. 

From a different point of view, one can imagine to keep A very small with 
respect to all the other energy scales present in the theory, making the new term 
a tiny perturbation of the system. Of course, nonperturbative effects can spoil 
this picture. 

The U{1) symmetry in standard BEG is related to the fact that the number 
of atoms is conserved. The breaking of this symmetry therefore is connected 
with the failure of this charge to be conserved. There are at least two possible 
physical implementations of this. 

A rather natural option is to have an open system. Goncretely, one could 
imagine to have a condensate which is able to exchange particles with some 
sort of reservoir, in such a way to preserve, on average, their number. Several 
settings in this sense could be conceived, e.g. with coupling with suitably tuned 
lasers. 

A second important case is the one in which the constituents themselves are 
some sort of collective degrees of freedom. This could be the case, for instance, of 
the excitations in the so-called quantum Heisenberg ferromagnet defined with 
a spin system. In this case, the fundamental operators give rise to effective 
degrees of freedom, called magnons, whose Hamiltonian, in general, is not U{1) 
invariant. For more details see [28], and references therein. 



4.1 The condensate wave- function 

We consider the dynamics generated by (|49|) . from which we want to extract 
the equation of motion for the condensate ip. The evolution of the mean field 
is easily determined in terms of the eigenstates \E) of the Hamiltonian H: 
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= -— V^V - ^i^p- \ip* + V + 2KnEiJ + nmsip*, (60) 
2m 

where m^; = (£'|x^|£^), — {E\x^x\E) encode the effect of the non-condensate 
atoms. This is the generahzation of the BogoUubov-de Gennes (BdG) equation 
for the condensate wave-function to the case A 7^ 0. 

In the standard case, if we have N particles in the condensate, the number 
density of the non-condensate fraction can be estimated to be of order 1/N 
with respect to the number density of the condensate, when the condensation 
mechanism is particularly effective. In this sense, the terms m, n are of order 
1/iV. This can be safely exported to our case, with the slight modification of 
the meaning of the number N, which does represent only the average number 
of particles in the condensate (see (|52|) ). 

At zeroth order in the 1/iV expansion, we have the generalization of the 
Gross-Pitaevski (GP) equation: 

d h? 

?/, = VV - mV' - A?/-* + kIt/-! V^- (61) 
The time independent homogeneous solution to the GP equation is 

nc=|^P = ^, (62) 



where we have fixed the phase of the condensate to be zero. In section | 
we will show that this is not an arbitrary choice, but rather a consequence of 
the situation we want to describe once we impose a stability condition for the 
quasi-particles. 

As in the standard case, we define the healing length ^ as the length scale at 
which the kinetic term is of the same order of magnitude of the self-interaction 
term in the Hamiltonian: 

Again, this length represents the spatial scale needed for the condensate to pass 
from the value ric = at the boundary of the region where it is confined to 
the bulk value ric (see [29]). Also in this case, the healing length represents the 
scale of the dynamical processes involving the deformation of the condensate 
wavefunction. This will have a crucial impact on the emergent gravitational 
dynamics in these systems. 



4.2 Quasi-particles 

The equation of motion for the particles out of the condensate is obtained by 
subtracting the equation for the condensate ((6T|l from the equation for given 
in (|49|) . We are interested in the propagating modes, so we neglect the self- 
interactions. We obtain: 

th^^X = -^^'X + (2«|^P - m)x + (^V-' - A)x^ (64) 
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Let us consider the ease of homogeneous condensate with density Uc given above. 
In this situation we have: 

d 

^^^X = -^V^X+(m + 2A)x + MX^- (65) 

If we decompose the field x in its plane wave components, we can rewrite this 
equation as 

d , _ h^e 

The mixing between a and due to the evolution in time becomes then appar- 
ent. We therefore pass to the quasi-particle operators ^(x) 



o h k 

ih—ak = -i^ak + (m + 2A)ak + lJ.a\. (66) 



^(x) = -^^Ske*-, (67) 

which are related to the particle operators through the Bogoliubov transforma- 
tion 

ak = a(fc)fek + /3(fc)fo-k, with a^{k)- (3"^ (k) = I. (68) 

The coefficients a, (3 are only functions of fc = |k|, since the condensate is homo- 
geneous and isotropic. The equation of evolution for the quasi-particles is then 
given by 

i/i^6k = £ik%, (69) 

with the energy 

The Bogoliubov coefficients are given by: 

where we have introduced the quantity 

A{k) = ^ + fx + 2X. (72) 

The high energy limit of these coefficients is: 

Hm a^{k) = 1, lim f3^{k) = 0, (73) 

k — ^oo k — >Qo 

which means that at large wave-number (and hence large momentum) , the quasi- 
particle operators coincide with the particle operators. This matches the behav- 
ior of the energy, which becomes just the energy of a non-relativistic particle 
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of mass m, just like a free atom. The dispersion relation ([70| suggests the 
introduction of the foUowing quantities: 



''sound 



(m + 2A)2 ^' > 



Here Csound plays the role of the speed of sound, while M plays the role of a rest 
mass for the quasi-particle. Since M is proportional to A, we clearly see that it 
is the term H\ that generates the mass of the quasi-particle. When A ^ 0, that 
is when H — > Hq, the quasi-particle becomes massless, i.e. a phonon, and the 
speed of sound reduces to the usual one in BEC. Perturbation theory, therefore, 
should be a viable strategy to compute the various physical properties of these 
systems. 

Notice that, in order to have a non- negative mass square term, and to avoid 
a tachyonic instability, we have to require A > 0. In standard BEC, one usually 
assumes that the chemical potential /i is positive: indeed if it were negative, 
there could not be any condensation. In our case, we can relax this requirement 
and obtain that /i > —A as a condition. In the following we consider /i > 0, in 
order to be able to consider the case in which the correction we are inserting 
is very small, without affecting dramatically the condensation. Indeed, it is 
easy to see that a condensation can take place even in a system with this soft 
U{1) breaking by checking the behavior of the two points correlation function 
G(x, y) = (>I'''^(x)>I'(y)). It is immediate to realize that, in the case of homo- 
geneous condensate, this correlation function describes long range correlations, 
since the mean field -0 is non- vanishing (cf. equation (1621) 1. 

M. is proportional to m, the mass of the atoms. Defining the ratio C — A//i, 
we introduce the function F{Q) 

= FicW = 4^^^m^. (75) 

Under our assumptions, we have that ( > 0. It is then straightforward to check 
that on this domain i^(C) is a monotonic (increasing) function and that 

F(0) = 0, lim F(C) = 1. (76) 

We conclude therefore that the mass of the quasi-particles A4 is always bounded 
by the mass of the atoms, M G [0, to). 

It is also interesting to notice that using the variable C, the speed of sound 

is: 

2 ^ l + KTlc , , 

For C small, we then have c^q„„^ « Kric/rn, which is the standard result, while, 
for C -> oo, c2^„„^ ^ 2KnJm. 
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4.3 The various regimes for the modified dispersion rela- 
tion 

Before moving on to the gravitational dynamics, let us discuss briefly the content 
of the dispersion relation ([70)1 for the quasi-particles, rewritten using Csound and 
M. 

( 

£{P) = + cloundP" + M^ctundj , (78) 

where we are using the obvious notation p ~ hk to simplify the shape of the 
expressions. Let us define the characteristic momenta pA, Pb and pc such that 

4 4 

Pa _ 2 2 Pb _ yv^2 4 2 2 „ yv^2 4 

/^Ul^ '-'soundP A-i if^rn^ sound! ^soundPC ' ^soundJ \'^J 

so that they are explicitly 

Pa = 4m2c2^„„^, p% = 2mMcl^^^a, p^ = A^^c^„„„^. (80) 
They are related through the relations 

pI ^ 2F{0pI ^ iF\C)p\. (81) 
If C ^ 1, which will be the regime we shall consider, we have also that 

PC <^PB<t:pA- (82) 
Taking into account (|82p . the characteristic momenta define different regimes: 

• li p ^ Pa, the term dominates, the dispersion relation (j78p is well 
approximated hy £ ^ p'^ /2m, we are in the trans- phononic regime. 

• If on the contrary we have pc ^ p ^ Pa, we can safely neglect the term 
of order p^, we are then in the relativistic regime since the dispersion 
relation ([78]) is weU approximated by f ~ {p'^^^lound + -^^^tound)^ ■ The 
quasi-particle is then relativistic, when the speed of sound Cs is playing 
the role of the speed of light. 

• If we are in the regime where p ^ pc, this means that the quasi-particle 
has a speed much smaller than Cgound, so that this is the Galilean limit of 
the relativistic regime. We are then dealing with a Galilean quasi-particle. 
The rest mass -Mcl^und provides the usual constant shift of the Galilean 
energy £ ~ Mc^^^^^ + p''/2M. 

4.4 The fluid description 

We have already seen that the standard Gross-Pitaevski (GP) equation describ- 
ing a BEC admits an interesting fluid interpretation, through the Madelung 
representation. We are considering now the GP equation given in (|6ip 

ot 2m 
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and we want to use the Madelung representation for the complex field i/}-. 

V = ^/^e-'^/^ (84) 

When replacing this into the GP equation, dividing by the phase and split- 
ting the resulting expression into the real and imaginary parts we obtain two 
equations: 

fic + V ■ (ricv) = --ricsmi—j , (85) 

= Kiuantum + - /i - A COS ^ - KTlc, (86) 

where we have introduced the velocity field v = —V9/h, and 

Kiuantum = ^ \/^' i^'^) 



is the familiar quantum potential term. These two equations, in the case A = 0, 
reduce to the usual form of the continuity equation and the Euler equation for a 
perfect fluid, once we neglect the quantum potential term. On the other hand, 
when A 7^ the i7(l) invariance is broken, and the number operator is no more 
conserved by the Hamiltonian evolution. 

It is interesting to see what happens when wc consider the case of homoge- 
neous condensates, d^ric = d^v'^ = 0. From the first equation we get: 

sin (^^^ =0<^e= y/i, I e Z. (88) 

This result implies that not only iJ is constant, but it actually vanishes. Inserting 
this result in the second equation wc obtain: 

^^^M + cosOtOA^ (89) 

K 

The analysis of the quasi-particle dynamics in the case of homogeneous conden- 
sates has shown that the case A < corresponds to a negative mass square term, 
i.e. tachyonic behavior: the energy of a quasi-particle would get an imaginary 
part leading to exponential growing and damping of modes. Since a choice of the 
phase of the condensate such as cos(?7r) = — 1 would be completely equivalent 
to a change of sign of A, thus leading to instabilities, it is clear that cos{Itt) = 1 
is required for the stability of the condensate. 



4.5 Gravitational dynamics 

The next step is the analysis of the inhomogeneous condensate, and hence the 
promised emergence of a gravitational dynamics. To simplify further the treat- 
ment it is better to consider the case of condensates which are nearly, but not 
exactly, homogeneous: this will correspond to the case of weak gravitational 
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field. This limitation is consistent with the formalism we are using. Indeed the 
mean field method certainly is not a good approximation in regions where there 
are large variations in density (see for instance vortex cores). 

In an asymptotically flat spacetime, in order to identify the Newtonian 
gravitational potential it is necessary to evaluate the non-relativistic limit of 
the geodesic equation in a weak gravitational field [3D]. In the asymptotic 
region there is a coordinate system such that the metric can be written as 
dtiu = Vfj.i' + h-fj.u- The tensor /i^j^ encodes the deviation from exact flatness, i.e. 
the gravitational field. In this regime, it is easy to show that the Newtonian 
gravitational field is identified with the component Hoq: 

$iv(x) = -i/ioo(x). (90) 

In the context of standard BEC (i.e. dealing with the non-linear equation ([T])), 
the quasi-particles travel in an emergent metric ds^ determined in terms of the 
homogenous condensate ip. 

ds^ = ^ [- (cs - v^) df - 2vidtdx' + Sijdx'dx^] , (91) 

where m is the mass of the atoms and and v depend on the properties of the 
condensate ip = y^e'^, through 

uric ^ 1 - 
Cs = , V — — vO. 

m m 

Considering that the condensate is homogenous, the density and velocity profiles 
become constant, i.e. respectively Uc — tioo, v — Voc- With the coordinate 
transformation, 

dT = dt, dX' = dx' ~ v'^dt, (92) 
the line element (|9ip is rewritten as: 

dsl^~cldT^ + dX^. (93) 

The condition of asymptotic fiatness for the acoustic spacetime can be trans- 
lated into the condition of asymptotic homogeneity for the condensate. We 
require then that only in a small region of space, in the bulk, the condensate 
deviates from perfect homogeneity. 

We consider therefore some small deviation from the asymptotic values of 
the velocity and of the density: 

JT-c = ^^00(1 + 2u(x)), v — Voo+w{x), with u ^ 1, W'^v. (94) 

This implies in particular a rescaling of the speed of sound. 

cl^^^cUl + 2u{x)). 
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The acoustic line element (lOTI) becomes then 



ds^ = J1^'J±^ (_(c2 _ y^)de - 2vd^'dt + J^dx^dx^) , (95) 
mcs noo 

where we have introduced a constant prefactor mCo^jn^ in order to have the 
conformal factor asymptotically normalized to one. Using (j94l) . together with 
the coordinate change ([9^ . the acoustic line element has the form: 

ds^ = dsi^ - ZuiXyl^dT'^ - 2w^{X)dTdX' + u{X)S,jdX'dX^ , (96) 

at first order in u, Wi. Consequently, we see that 

Ko[X) - ^ic^^u{X), (97) 

so that the gravitational field is encoded in the number density perturbation of 
the condensate wave-function tp^ 

^n{X) = \cl^u{X), (98) 

while it is independent from velocity perturbations, which therefore can be 
discarded. 

This result allows a simplification in the choice of the physical situation: 
it is enough to discuss the case in which the condensate wavefunction has a 
constant phase, while its modulus slightly deviates from perfect homogeneity. 
It is convenient to introduce the parametrization: 

\ \ 1/2 



where u(x) is a dimensionless function and it is assumed to be very small. In 
practice, we will assume that it is associated with a localized inhomogencity of 
the condensate. At infinity we ask that u 0. Notice that the wavefunction 
is real, due to the fact that we are really interested in number density 
fluctuations, and not on fluctuations in the velocity profile. This simplification 
reduces the number of independent functions without making the system trivial, 
as it will be shown. 



4.6 The gravitational potential for the quasi-particles 

Having discussed the setup in which we are working, we can reconsider the 
quasi-particle dynamics in this new system. It is not necessary to recover some 
sort of acoustic metric. While it is certainly interesting, given that the quasi- 
particle are massive and that we are interested in the nonrelativistic limit, the 
notion of acoustic metric is of little interest, in this particular case: it has been 
used just as a guide to isolate a candidate for the gravitational potential. 

The first step, then, is to see if there is a term in the equation of motion for 
quasi-particles (|64p . which can be identified as an external potential term. This 
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will allow us to check the conjecture that it will be given by the number density 
perturbation u(x), as well as the precise coefficient relating it to the familiar 
Newtonian potential (having dimensions of the square of a velocity) . The next 
step will be to take this potential and plug it into the Bogoliubov~de Gennes 
equation (j60p . which is describing its dynamics. 

To identify the Newtonian potential, the diagonalization of the Hamiltonian 
in ((64|) for the field x must be done again, including now the fluctuations of the 
condensate wavefunction. In this case, the diagonalization procedure is more 
involved: we have to deal with the non-commuting operators and u. We 
can not perform it in an exact way. However, we are interested in the Galilean 
regime for the quasi-particle spectrum, when pc ^ P (cf section . It is then 
a reasonable approximation to neglect all the terms involving the commutators 
[p^/2to, m(x)], which are largely suppressed (with respect to the other terms 
appearing in the equations) by the mass of the atoms and from the smallness 
of u(x). 

With these simplifying assumptions, the Hamiltonian for the quasi-particles 
in the non-relativistic limit is 



H„ « McLw + 2 (^i±^l(^i±i^.(x), (100) 

SOU7ld 



where the mass of the quasi-particle M and for the speed of sound Cs are given 
in ([7^1) . We first recognize the constant shift Alc^^^^^ of the energy due to 
the rest mass in the Galilean regime. This term is not affecting the discussion 
in any way and can be subtracted without physical consequences. The term 
proportional to ■ti(x) can be clearly interpreted as an external potential. If we 
want to identify it with the gravitational potential $grav, we need to have 

MC^ ^ = A^^grav ^ $grav(x) = ^ ^ ^w(x), (101) 

where A4 is the mass of the quasi-particles. Note that this identification is 
formal, and relies on the way in which the gravitational potential enters the 
Schroedinger equation for a non-relativistic quantum particle. We should always 
work with u: our definition of $giav is dictated from the analogy we want to 
make with Newtonian gravity. For instance, we see that this definition becomes 
singular when we deal with massless quasi-particles, i.e. when A ^ 0. This must 
be expected: when A vanishes the quasi-particles become massless phonons, for 
which the coupling to a Newtonian gravitational potential cannot be defined in 
terms of their mass density. 



4.7 The modified Poisson equation 

Now that we have identified a candidate for the Newton potential 3'grav from 
the quasi-particles dynamics, we need to check that it satisfies some sort of 
Poisson equation. Since the gravitational potential is deduced from ^ - as 
small deviations from perfect homogeneity (cf. ([99]) ) - the Poisson equation 
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should be deduced from the BdG equation ([SU]) . With the natural assumption 
that the potential is reacting instantaneously to the change of distribution of 
matter, we can neglect the time derivative and (|60p becomes 



- 2(^I + A) m(x) = 2k n(x) + -m(x) 

2m / V 2 



(102) 



We have seen in section ITT] that the terms m(x) and n(x) are functions of the 
atoms X outside the condensate and therefore of the quasi-particle </), through the 
Bogoliubov transformation ([55)1 . Therefore they can be interpreted as the source 
in the (modified) Poisson equation. We examined different types of source: 
either localized particles or plane-waves. We shall discuss here only the first 
kind of sources and forward the reader to |3T] for the second kind. 



4.7.1 Localized sources 

The most natural source to consider for the Poisson equation is a single quasi- 
particle at a given position xq. However, point-like distributions give rise to 
divergencies. We consider therefore a quasi-particle which is localized around 
the point xq, with a non-zero spread to regularize these divergencies. We con- 
sider a quasi-particle in a state of the form: 



ICxo)- / dVxo(x)^t(x)|l]), 



with 



Cxo encodes the spreading of the particle around xq since 

(Cxol0nx)0(x)|Cxo) = ICx„(x)p. 



d3x|Cxo(x)|2 = l^(C.JCxo) = l. 

(103) 



(104) 



We can now determine the value for the anomalous mass m and anomalous 
density n when the quasi-particle is in the state |Cxo)- explicit calculation 
gives 



n(x) 



dV(x-z)Cxo(z) 



d3z5(x- z)Cxo(z) 



l5]/?2(k),(lG5) 



m X 



2(^1 d3zig(x-zi)Q(zi)^ (^J d3zi/(x-Z2)Cx„(z2)) +-^E"(k)/3(k)(106) 



where we have introduced the functions /, g depending on the Bogoliubov co- 
efficients a and f3 



/ W = ^ E «(k)e*■^ 5(x) 

k k 

The quantities and mo with 

"f^ = -^E/5'W, mn = -^5]a(k)/3(k). 



V 



(107) 



(108) 
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are vacuum contributions independent from the presence of actual quasi-particles. 
They are related to the inequivalence of the particle and quasi-particle vacua, 
and it can be easily seen that: 

The functions /, g encode the fact that quasi-particles are collective degrees of 
freedom and therefore intrinsically some non-local objects. This non-locality 
is due to the Bogoliubov transformation ((68|) . Quasi-particles and atoms (i.e. 
local particles) coincide only if we have a(k) = l,/3(k) = 0, and therefore 
/(x) = (5'^(x), while ^(x) = 0. Since this is not the case, the anomalous mass 
and the anomalous density will show an intrinsic non-locality. The spreading 
characterized by |Cxo) encodes some extra non-local effect, introduced by hand 
for regularization purposes. Therefore this feature is not as fundamental as the 
non-locality introduced by the Bogoliubov transformation. 
The equation (|102p becomes then: 

- 2(/i + A)^ u{x) = 2k (^n{x) + ^m(x)^ + 2k {nn + ^mo^ , (110) 

where we have introduced the quantities 

n(x) = n(x) - nn, m(x) = m(x) - mo, (HI) 

which represent the contribution of actual quasi-particles to the anomalous den- 
sity and anomalous mass, respectively. By dimensional analysis, the terms n, m 
have the dimensions of number densities. Since in Newtonian gravity the source 
for the gravitational field is a mass density, we introduce the mass density dis- 
tribution: 

1 , 



Pmatter 

) = M l^n(x) + 2"^W j • (112) 
With this definition, we can rewrite (jllOp as an equation for the field <i>grav: 

- ^ $g,a^ = 47rGp,„attcr + A, (113) 

where we have defined 

^ ^(m + 4A)(m + 2A)^ 2^(A^+4A)(Ai + 2A) , 1 
^ = 4.n2^A3/2(^ + A)V2' ^ = WX ("^^ + 2"^^^^') 

Am{^. + A) ^ ' 

This particular choice of notation is motivated by the comparison of (|113p with 
the Newtonian limit of Einstein equations with a cosmological constant. 

For this reason, we can identify these three quantities as the analogous of 
the Newton constant, the analogous of the cosmological constant and a length 
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scale which represents the range of the interaction, as we are going to discuss 
below. 

To get a better grasp of the physics of the modified Poisson equation (|113p . 
we can look at its solution for a given distribution of quasi-particle Pmattor- 
As it is well known, a solution for the equation 

W^-^^^x)^iTTGMS^{x-z), (116) 

is given by the Yukawa potential 

4>v(x;z) = ^--^. (117) 

On the other hand, a solution for the equation 

1 



- j $(x) = A, (118) 

is just given by the constant solution 

$A = -L^A. (119) 

Notice the peculiarity of this solution. It does not give rise to a gravitational 
acceleration since the gradient is trivially zero. Therefore the only effect of this 
term is to shift the overall density of the condensate. 

The linearity of equation pi3p allows us to use these results to write down 
a solution for a generic distribution of matter (i.e. quasi-particles) as 

*grav(x) = yPmattcr(z)*Y(x;z)d3z + $A. (120) 

Solutions of (|113p are therefore constructed from the Yukawa potential smeared 
out due to the non-locality of the quasi-particle (with an extra global shift due 
to the cosmological constant). The Yukawa potential is typically encoding some 
short range interaction, characterized by the scale L which is simply related to 
the heahng length ([63]) . 

= ^. (121) 

Although this a very short range for the gravitational interaction, this outcome 
should not come as a surprise. In fact, the healing length (c.f. (|63p) characterizes 
the typical length over which a condensate can adjust to density gradients. Since 
density inhomogeneities encode the gravitational interaction, one should expect 
them to be damped over a distance of the order the healing length. 

In the context of relativistic field theory, the short interaction scale for grav- 
ity would be translated in a massive graviton, with mass given by 



sound ^ ' 
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We can then compare the masses of the quasi-particles , graviton Mgrav and 
atoms TO, 

< M < m < V2m < Mgrav < '2m, (122) 

which shows the hierarchy of the energy scales present in this system. We 
notice that the graviton is then always more massive than the quasi-particles, 
and that this interaction is of very short range, since the ^ is much shorter than 
the acoustic Compton lengtli3 of the quasi-particles. In particular, we cannot 
tune the parameters of the system in such a way to make Mgrav arbitrarily 
small, in order to be closer to reality. 

4.8 Lessons 

In an analogue gravity model based on a BEC system, the degrees of free- 
dom are separated into the atoms that condense and the ones which do not. 
Quasi-particles are then collective degrees of freedom constructed from the non- 
condensed atoms. The dynamics of the quasi-particles can be described, in a 
given regime, in terms of the propagation of particles over an effective curved 
spacetime metric, which is a function of the density Uc and the velocity profile 
V of the condensate. In this sense, it is natural to expect that gravitational 
degrees of freedom are encoded in the condensate. Dynamics of the latter is 
encoded in the BdG equation ([60|) . which is essentially Galilean. Hence, we can 
not expect to recover the Einstein equations in this context [l^- Nevertheless, 
one can still try to interpret ([60)1 as some sort of Poisson equation for some type 
of Newtonian gravity. 

However, quasi-particles are massless in usual BEC systems and hence they 
cannot be considered as sources for the gravitational field in the Poisson equa- 
tion. We introduced therefore a new term Hx in the dynamics of the BEC which 
softly breaks the U{1) symmetry and consequently, as we showed in section l472l 
generates a mass gap for the quasi-particles. We showed explicitly that the 
presence of this small symmetry breaking term does not prevent a condensation 
from happening and still allows a mean field description (which is sufficiently 
accurate for our purposes). Then, following the usual general relativistic argu- 
ment, we have argued, in section (|4.5p . that the Newtonian potential has 
to be related to small inhomogeneities in the condensate density (while pertur- 
bations in the velocity profile do not contribute at first order as gravitational 
degrees of freedom). This conjecture, based on the analysis of a standard BEC 
system, was then confirmed by a specific analysis of the modified BEC dynamics 
for an almost homogenous condensate. 

The end point of this investigation can be then summarized in the following 
two equations 

i? = Xa= -7WV$grav, (123) 

- -^^ *grav = ^ttGnP + A, (124) 

^We are using Cgaund instead of cngi^t to define all these scales. We have to use the natural 
units for a hypothetical phononic observer. 
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where M is the mass of the quasi-particle acquired via the soft U{\) symmetry 
breaking induced by ((48|) . L is proportional to the heahng length, A plays the 
role of the cosmological constant and Gn is an effective coupling constant that 
depends on the condensate microphysics and the form of the matter source. 

For what regards the latter we have considered two cases: a localized quasi- 
particle state and a set of plane waves. In the first case the analogue Newton 
constant is indeed momentum and position independent and the solution of the 
modified Poisson equation p24p , has the form of a smeared Yukawa potential. 
The smearing is due to the fact that quasi-particles are intrinsically non-local 
objects, being collective degrees of freedom. When considering plane-waves 
as sources, we have instead that, due to the momentum dependence of the 
Bogoliubov transformation. Gat is running with the momentum and the solution 
for the gravitational potential is a constant (albeit a different one for different 
momenta). One should however be careful: while it is common in quantum field 
theory (QFT) to encounter the notion of running coupling constants, the origin 
of the running here is rather peculiar. Indeed, in QFT the running is due to 
quantum corrections to the tree level/classical action, here the running is due 
to the inequivalence between the ground state of the Fock spaces of atoms and 
quasi-particles. Paraphrasing what has been done in the context of emergent 
geometry, where the notion of "rainbow geometry" has been introduced, we 
could speak about "rainbow dynamics" . 

We have also obtained naturally a cosmological constant in the model: vac- 
uum gravitates, even though in a very peculiar way. It is induced by the terms 
(r2|x^xl^)) (^lxxl^)i where is the state with no quasi-particles. It is entirely 
due to the (unavoidable) inequivalence between the quasi-particle vacuum and 
the particle vacuum and cannot be put to zero just tuning the parameters. It 
represents an interesting alternative to known mechanisms to generate a cosmo- 
logical constant (see also [13] for similar ideas about the nature of the vacuum 
energy in condensed matter systems). 

Let us compare this point with the standard cosmological constant problem. 
If one thinks to the cosmological constant is generated by the zero point energy 
associated to each mode, in the case of BECs a naive expectation would be 
that the cosmological constant term would be set at the characteristic UV scale 
represented by the inverse of the healing length. This would lead to a very large 
vacuum energy. However, we have seen that the cosmological constant, in this 
specific case, is linked to the so called depletion factor (i.e. the ratio between the 
non-condensed fraction and the condensed one), and hence naturally suppressed 
by a factor \/N in the expansion in the inverse of the number of atoms. It is 
the condensation mechanism itself guaranteeing a naturally small cosmological 
constant. Of course, the worse is the condensation, the larger is the number 
of atoms out of the condensate and hence the bigger will be the cosmological 
constant. This fact has a rather peculiar interpretation. 

To have an emergent Lorentzian structure, it is essential that there is a mean 
field which represents the emergent spacetime. This was encoded in the splitting 
of the field operators ^ « V + X- However, for a mean field approach to make 
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sense, fluctuations around the mean field should be relatively small: 



< 1. (125) 



Therefore, the smallness of the analogue of the cosmological constant term in 
the BEC is deeply intertwined with the very definition of mean field, i.e. how 
good is the picture of quasi-particles moving in a classical Lorentzian spacetime. 
Of course this poses the question of whether this mechanism do apply also for 
spacetime in which we actually live. 

In conclusion, BEC as an analogue model for gravity presents many differ- 
ences with a realistic gravity theory as we expected. We do not get general 
relativity in a condensate. However 

• there is an emergent Lorentzian metric describing the propagation of the 
phonons; 

• there is a mismatch between microlocality and macrolocality due to the 
nonlocal nature of the phonons; 

• there is an emergent Newtonian gravitational theory, which is very short 
range; 

• in this theory vacuum gravitates; 

• the source term for the gravitational field inherits the nonlocality proper- 
ties of the phonons; 

• the cosmological constant is naturally small provided that the depletion 
factor is small. 

Despite the limited applicability of the results, the suggestions for realistic 
theories of gravity, in particular for quantum gravity and the role of locality, and, 
perhaps more interestingly, for the cosmological constant problem are definitely 
worth of further investigation. 



5 Emergent gravity: the role of symmetries 

The discussion of the BEC model has shown that by limiting the analysis to 
condensed matter systems there are rather strong constraints on the kind of 
gravitational models which is reasonable to simulate. For example, in the above 
investigation we had a single scalar field: it would be interesting is to see what 
happens if several different species are present. In that case, besides the issue 
of having a short range rather than a long range interaction, also the coupling 
to the gravitational field must be carefully discussed. Indeed, in order to have 
some sort of equivalence principle, all the fields must be coupled to the grav- 
itational field in the same way. The natural setup to discuss these issues is 
the 2-BEC model [311 133] ■ in fact in this case one could treat a multi-particle 
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system whose richness could allow a closer mimicking of Newtonian gravity 
with a long range potential. However, the fact that emergent gravity has to be 
Newtonian in a BEC-based analogue model seems to be unavoidable since the 
gravitational potential depends on the condensate, which is typically described 
by non-relativistic equations. A possible way to avoid this issue is either to 
consider relativistic BEC [Ml [S^ (however in this case we would still expect to 
get only some type of scalar gravity), or to change completely paradigm and 
identify gravity not as the condensate but as linked, together with matter, to 
the perturbations around the condensate. We will consider later this second 
point of view in a different model. 

Furthermore, there is another important issue that needs attention. In our 
treatment we neglected the quantum potential, i.e. we have deliberately worked 
in the hydrodynamic limit of the theory, carefully avoiding the issue of the 
breakdown of acoustic Lorentz invariance in the system at suitably high energies 
of the quasi-particles. Presumably the breakdown of this emergent spacetime 
symmetry, namely local Lorentz invariance, will be linked also to some relevant 
regime change in the gravitational dynamics (which is anyway affected by the 
presence of a Lorentz symmetry breaking scale, the healing length, which ends 
up setting the graviton mass scale) . Should we take Lorentz symmetry breaking 
as a crucial ingredient of the emerging gravity paradigm or as an accident of 
the condensed matter analogue models? In the first case, how the breakdown of 
such spacetime symmetry affects the symmetries of the gravitational dynamics 
and in particular diffeomorphism invariance? Furthermore, does this imply that 
an emergent gravity scenario should give up the relativity principle and bring 
us back to Newton's absolute space and time? 

In order to explore these issues we can start investigating the role of Lorentz 
invariance in emergent gravity scenarios by considering the most well known 
"no-go theorem" against them, i.e. the so called Weinberg- Witten theorem 

m- 

5.1 LIV and emergent gravity: the Weinberg— Witten the- 
orem 

The idea of having the graviton as a composite particle/emergent field is cer- 
tainly a fascinating idea. However, there are limitations to what it is possible to 
do. In particular, there is a theorem, due to Weinberg and Witten [36], which is 
often presented as a crucial (fatal, in fact) obstruction for a successful emergent 
gravity program. 

The theorem states precise limits for the existence of consistent theories with 
massless particles. It has two parts, and it says that (quoting from [36] ): 

1. A theory that allows the construction of a Lorentz-covariantly 
conserved four-vector current J'' cannot contain massless par- 
ticles of spin j > 1/2 with nonvanishing values of the conserved 
charge / J^Sx. 
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2. A theory that allows for the construction of a conserved Lorentz 
covariant energy-momentum tensor 6^^'^ for which J O^'^dPx is 
the energy-momentum four- vector cannot contain massless par- 
ticles of spin j > 1. 

For a careful discussion of the proof of the theorem, and for references, see 
[37]. For additional comments, see [551 139j . 

Crucial ingredients for the proof of this theorem are Lorentz invariance and 
the nonvanishing of the charges obtained from Lorentz covariant vectors and 
tensors. Interestingly, the gauge bosons like the gluons and the graviton are not 
forbidden since the current for the gluons is not Lorentz-covariant conserved, 
and the graviton does not possess a covariant stress-energy tensor (but rather 
a pseudo-tensor). 

This theorem, then, poses rather strong constraints on the possible theories 
that can be built in Minkowski spacetime. Of course, gravity is not just the 
theory of a spin-2 particle in Minkowski spacetime. Nevertheless, it surely makes 
sense to consider the linearized theory in sufhciently small neighborhoods. In 
this limit, then, the theorem does apply. 

With this caveat in mind, we can say that in an emergent gravity pro- 
gram this theorem must be taken appropriately into account and appropriately 
evaded. There are (at least) two "obvious" way out: 

• allow for Lorentz symmetry breaking, or 

• make the spacetime manifold to emerge as well. 

The first option is rather straightforward, and it is essentially what could be 
pursued within scenarios like the one considered in analogue models, in which a 
preferred time function is specified. However, there is apparently also a (concep- 
tually high) price to pay: a step back from Minkowski spacetime to the notions 
of absolute space and time. Moreover, and most importantly, there is the issue 
of recovering a low energy approximate Lorentz invariance. We shall come back 
later on these issues. 

The second option is probably the most viable, conceptually appealing, but 
most demanding in terms of new concepts to be introduced. If no reference is 
made to a background Minkowski spacetime, but rather the graviton emerges in 
the same limit in which the manifold emerges, then there is no obvious conflict 
with the Weinberg- Witten theorem. Simply, what is called the gauge symme- 
try in terms of fields living of spacetime is the manifestation of an underlying 
symmetry acting on the fundamental degrees of freedom in the limit when they 
are reorganized in terms of a spacetime manifold and fields (gauge fields and 
gravitons in particular). There are already two examples of this possibility, 
namely matrix models and quantum graphity models. In both cases, the very 
notion of spacetime manifold is immaterial for the foundations of the theory. 
The manifold and the metric are derived concepts, obtained in precise dynam- 
ical regimes of the theory. The interested reader can find additional comments 
and references in [TUl [JD] . 
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The bottom line of this very concise overview of the WW theorem is clear: 
to obtain a realistic model of emergent gravity one must ask for very special 
mechanisms to be at work in the model. Without these, the theory would not 
be able to give a meaningful limit. 

5.2 Why breaking Lorentz symmetry might be good 

Specific hints of LV arose from various approaches to Quantum Gravity. Ex- 
amples include string theory tensor VEVs [JT], spacetime foam semi- 
classical spin-network calculations in Loop QG [43! , non-commutative geom- 
etry [HI H51 |3nj , some brane- world backgrounds i47j and condensed matter ana- 
logues of "emergent gravity" [T^lll^ • Although none of these calculations proves 
that Lorentz symmetry breaking is a necessary feature of Planck scale physics, 
they did stimulate research aimed at understanding the possible measurable 
consequences of LV [49l [50l [51]. Furthermore, recent investigations strongly 
suggested that an high energy breakdown of Lorentz invariance might strongly 
improve the renormalizability of field theories [521 153| including gravitation [54] . 

This seems in close analogy with what we naively expect from analogue 
models scenarios like the BEC one. There, in fact, the renormalizability of 
the fundamental field theory, a non-relativistic A** action, is preserved by the 
Bogoliubov transformation leading to the quantum field theory of the quasi- 
particles. Should we take these results as a strong hint that Lorentz symmetry 
breaking should be a part of any working emergent gravity scenario? It is 
at this stage unclear if we can be that bold. Surely one open issue is the 
naturalness of theories endowed with Lorentz symmetry breaking. In general, 
radiative correction, or alternatively renormalization group running, lead to a 
dangerous "percolation" in the infrared regimes of the Lorentz breaking [55|, 
[56t [57] , something strongly constrained by current observations [49l (50] [5T] . It 
seems that in order to solve this issue some sort of custodial symmetry would 
be needed (see e.g. the related discussion in [5B]) but no satisfactory solution 
up to date has been found. 

In addition to the above, more technical, issue it is also clear that a more 
philosophical one is also present. As we said, many find quite unappealing the 
idea to give up the relativity principle and go back to a preferred system of 
reference. In this sense, however, some comments are in order on the relation 
between Lorentz invariance and relativity. 

5.3 Relativity beyond Lorentz? 

Lorentz invariance of physical laws relies on only few assumptions: the prin- 
ciple of relativity, stating the equivalence of physical laws for non-accelerated 
observers, isotropy (no preferred direction) and homogeneity (no preferred loca- 
tion) of space-time, and a notion of precausality, requiring that the time ordering 
of co-local events in one reference frame be preserved [SSI IMl (SO] ■ In this sense 
a breakdown of Lorentz invariance does not necessarily imply a breakdown of 
the relativity principle. For this reason, it is worth exploring an alternative 
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possibility that keeps the relativity principle but that relaxes one or more of 
the above postulates. Such a possibility can lead to the so-called very special 
relativity framework [61] , which was discovered to correspond to the break down 
of isotropy and to be described by a Finslerian-type geometry [521 IMl IM] ■ In 
this example, however, the generators of the new relativity group number fewer 
than the usual ten associated with Poincare invariance. Specifically, there is an 
explicit breaking of the 0(3) group associated with rotational invariance. 

One may wonder whether there exist alternative relativity groups with the 
same number of generators as special relativity. Currently, we know of no such 
generalization in (commutative) coordinate space. However, it has been sug- 
gested that, in non-commutative spacetime, such a generalization is possible, 
and it was termed "doubly" or "deformed" (to stress the fact that it still has 10 
generators) special relativity, DSR ^5]. Unfortunately, the various DSR can- 
didates face in general major problems regarding their physical interpretation 
(e.g. the so called "soccer ball" problem [65]). 

Finally, it is a logical, and rather simple, possibility that a Lorentz symmetry 
breakdown could be signaling an interpolation from a relativity group to another 
one, for example two special relativity groups characterized by different limit 
speeds or between a Lorentzian and an Euclidean Poincare group. This second 
possibility is quite appealing because it would allow to consistently introduce a 
minimum length without reducing the generators of the Poincare group as well 
as it would give a natural meaning to the breakdown of Lorentz invariance as 
we know it by associating it with the emergence of time. 

In what follows we shall pursue this route and try to built up a model of 
emergent gravitational dynamics where time and gravity will emerge from a 
Euclidean manifold endowed with a set of fields. 

6 Emerging time and scalar gravity 

The discussion of the BEC model has shown that by limiting the analysis to 
condensed matter systems there are rather strong constraints on the kind of 
gravitational models is reasonable to simulate. Therefore, we will leave the 
arena of analogue models and we will present a toy model in which a number of 
different issues can be addressed. In particular, we will focus on two of them. 

First of all, given that analogue models are condensed matter systems, the 
notion of time is inherited from the time of the laboratory. This structure 
percolates from the fundamental level of Galilean spacetime where the atomic 
dynamics takes place onto the effective dynamics of the propagating degrees of 
freedom, e.g. phonons. 

There is an exception to this: by a careful tuning of the coupling constants, 
it is possible to make the scattering length of a BEC a negative quantity, and 
hence produce an effective dynamics for phonons which is Euclidean and not 
Lorentzian. Of course, this is a very interesting phenomenon which gives the 
possibility of studying a rather exotic class of phenomena related to signature 
change events |66[ 167] . Nonetheless, this phenomenon is rather simple, from a 
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conceptual point of view: while time is present at the fundamental level as a 
definite structure (i.e. a preferred class of foliation of the spacetime manifold), 
this latter is hidden for the emergent system. 

It would be much more interesting to understand how it is possible to do 
the opposite, i.e. whether it is possible to emerge time out of a timeless system. 
We will show a simple toy model in which this happens [68j . 

The second point concerns the issue of diffeomorphism invariance. Diffeo- 
morphism invariance is one of the distinctive features of General Relativity. In 
an emergent gravity program, it is crucial to understand how it will be possible 
to emerge it as well. By elaborating further on the toy model used to emerge 
time, we will show how a diffeo-invariant theory, namely Nordstrom theory for 
scalar gravity, can be extracted from this particular example. 



6.1 Emergence of time 

Let us assume that a fundamental unknown theory gives rise in some large 
number limit to simple structures such as R** equipped with the Euclidean metric 
5^"^, and a set of scalar fields ^'^(a;^), i — 1, ...,N {x^ £ M**) with their Euclidean 
Lagrangian C. Since we do not know this fundamental theory, we choose such 
Lagrangian to be of the simple shape0 



(126) 



with Xi = 5'^'^ d^'^ idv'^ i. It is easy to see that this Lagrangian is invariant 
under the Euclidean group ISO{4). The equations of motion are then simply 
for a given field ^'i 



dF 
dXl 







/ d^F 



dF 



(127) 



Let us now consider a specific solution of the above equations of motion, 
and perturbations (pi around it. For — ipi + (pi, the kinetic term Xi becomes 
then 

Xi + 5Xi, with Xi = S^" dfj,ipid^ipi and 

SX,^2df,yj,df'^, + df,ip,d^'^,. (128) 

We intend now to identify some specific F such that the Lagrangian for the 
perturbations (pi is invariant under the Poincare group ISO{S, 1). To determine 
the Lagrangian for the perturbations (pi, we expand p26p using (|128p . 



F{Xi, .., Xn) 



F{Xi,..,Xn) + 



d^F 



^ dX, 



2 ^ dXjdXk 

jk J 



5X, 
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_^X,5Xu , Q^Qx,dXM 

jkl ■> 



8Xj,8Xk6Ximi) 



■^We could also consider a dependence on crossed terms of the kind h>^" d^^^ idv^ however 
this is not changing the final result. 
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The first term F{Xi, ..,Xn) is the Lagrangian for the classical solution ipi. The 
second term, the one linear in SXj, contains a term linear in d^tpi, which is zero 
on shell. We can also identify the quadratic contribution for dfj^ipkdyipk'- 



X 



for k = l, df^ipkd^ipk 



OF 



^dXk 



5^" + 



1 d'^F 



X 



2 [dXkf 



(130) 
5^Vfc9''7/.fe ) .(131) 



X 



The contribution (|130p introduces some mixing between fields in the kinetic 
term. To simplify the analysis, we demand that they cancel, which puts a 

= 0, if fc ^ Z. A specific solution is 



dXkdXi 



constraint on the choice of F, i.e 
then 

F{X^,..,XN)^h{Xi) + ... + fN{XN). (132) 

We can identify in (113111 the effective or emergent metrical for each field (pfe, 
(taking into account ()132p ) 



9k 



dfk 



dXk 



Xt 



cPfk 



[dXkf 



(133) 



Xi, 



Since a priori fi ^ fj and ipi ^ if * ^ J j we are dealing with a multi-metric 
structure: each field sees its own metric. However, we can enforce a mono-metric 
structure by constraining the solution ■0^ and the derivatives of fk at Xk to be 
independent of k 

fk^f, ^k^^, Vfc. (134) 

So far we have just shown that the perturbations around a solution of the 
field equations on a Riemannian manifold can propagate, for suitably chosen 
Lagrangians, on an effective geometry which is not the fundamental one, 
but rather a rank 2 tensor constructed from it and partial derivatives of the 
chosen background solution. Note that, in order for this to be possible, it was 
crucial to have a starting Lagrangian with non-canonical kinetic terms as it can 
be clearly evinced by the second contribution to the metrics in equation ([133]). 
As a next step, we show now how for some solutions of the equations of motion, 
such effective metric can be of pseudo-Riemannian form. In fact, we can even 
ask that the metric (jl33p is the Minkowski metric rj^^. This will put some 
constraints on the derivative of /, evaluated sA X = d^ipdf^ip. 

In order to do so, we shall need to specify a particular solution, ^5, of the 
equations of motion. Let us take it to be an affinc function of the coordinates, 
= a^Xfi + f3. It is easy to check that this is indeed a solution of our field 
equations p27p . Moreover, thanks to the S'0(4) symmetry, we can always make 
a rotation such that 

■iP = axo + l3. (135) 



•^Actually, we show here the inverse metrics from which the actual metrics can be derived 
once invertibiUty conditions are imposed. In our case of interest, this will always be true. 
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The choice of the coordinate xq is completely arbitrary, what only matters is 
that there is one coordinate which is pinpointed. Finally, we ask for the metric 
to have the signature (— , +, +, +)■ This puts some constraint on the value of 
the derivatives of / 



dX 

dl_ 

dX 



1 dy 

y+2 {dXr 

. 1 dy 



< 0, 



X 2 {dxy 



d''ij}d°"ip > a = 1,2, 3 
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which using (I135|) imply 
dX 



X 



a 
T 



[dxy 
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dX 



> 0. 



(136) 



(137) 



X 



Note that, due to the choice of a solution of the form (|135p . the conditions p36p 
are not only implying a pseudo-Riemannian signature but also the constancy 
of the metric components, which hence can be easily rescaled so to take the 
familiar Minkowskian form diag(— 1, +1, +1, +1). 

Of course, there are many possible choices of f{X) and a which can fulfill 
the above requirements. For example, we can pick up the specific combination 



f{X) 



-x' + x. 



1 



1 



3<" <2 



(138) 



However, in what follows we should not make use of any particular form oi f{X) 
and a and simply assume that they are such that (|137p are satisfied. 

To summarize, since g'^'^ = rj^'^ ^ Wk, the (free) perturbations (pi are propa- 
gating on a Minkowski space, even though the fundamental theory is Euclidean 
(c.f. (|126p ). At this point few remarks are in order. 

So far, our theory does not posses any fundamental speed scale. This is 
natural since the fundamental theory is Euclidean. At this level, there is no 
coordinate with time dimension and therefore one cannot define a constant 
with speed dimension. The invariant speed c, which will relate the length xq 
to an actual time parameter t, could be determined experimentally by first 
introducing a coordinate with time dimension (as it would be natural to do 
given the hyperbolic form of the equations of motion for the perturbations) and 
then by defining c as the signal speed associated to light cones in the effective 
spacetime0. 

Second, a comment is due about our choice of the background solution 
around which we have considered the dynamics for perturbations. It is ob- 
vious that within our model this choice is arbitrary. It simply shows that there 
are some background solutions ip for which a pseudo-Riemannian metric can 



* Noticeably, a similar situation is encountered in the von Ignatowsky derivation of Special 
Relativity 1581 where, given a list of simple axioms, one derives the existence of a universal 
speed, observer independent, which is not fixed a priori to be the speed of light but has to be 
identified via actual experiments. 
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emerge. Obviously, diflFerent background solutions could lead to alternative 
metrics, e.g. one could also obtain the Euclidean metric S^^, (for example if if) is 
constant), a degenerate metric or more complicated structures according to the 
possible solutions f/'. While it is conceivable that in a more complicate model we 
could have some mechanism for selecting the specific background solution that 
leads to an emergent Lorentzian signature, it is not obvious at all that such a 
feature should be built in the emergent theory. In fact, one generally minimizes 
an energy functional to select the ground state of the theory. However, when 
looking at Lorentzian signature emergence starting from an Euclidean set up as 
in our model, there is no initial notion of time and hence no energy functional 
to minimize. It is therefore unclear how a ground state could be selected from 
within the emergent system. 

On the other hand, it is also conceivable that the actual background solution 
in which the initial system of fields (|126p emerges from the fundamental (pre- 
manifold) theory, can be depending on the conditions for which the "conden- 
sation" of the fundamental objects takes place. In this sense, the right ground 
state or background solution would be selected from minimizing some functional 
defined at the level of the atoms of space-time. To use an analogy, the same 
fundamental constituents, e.g. carbon atoms, can form very different materials, 
diamond or graphite, depending on the external conditions during the process 
of formation. Similarly, in a Bose-Einstein condensation the characteristics of 
the background solution (the classical wave function of the condensate), such 
as density and phase, are determined by physical elements (like the shape of 
the EM trap or the number and kind of atoms involved) which pre-exist the 
formation of the condensate. 

In conclusion, we have identified the fundamental Lagrangian so that the 
perturbations Lpi have a kinetic term determined by the Minkowski metric. 

£cff('^i, --Vn) = ^rfdf^ipidi.ipi. (139) 

i 

In this sense, we have a toy-model for the emergence of the Poincare symmetries. 
This construction can be seen as a generalization of the typical situation in ana- 
logue models of gravity [12] where one has Poincare symmetries emerging from 
fundamental Galilean symmetries [T^. However, let us stress that in our case no 
preferred system of reference is present in the underling field theory given that 
the fundamental Lagrangian is endowed with a full Euclidean group IS 0(4:). 
Moreover, the emergence of a pseudo-Riemannian metric is in our model free of 
the usual problems encountered in the context of continuous signature change 
(e.g. degenerate metrics) given that the former arises as a feature of the dynam- 
ics of perturbations around some solution of the equations of motion. Similarly 
one can see that the invariance under Lorentz transformations is only an ap- 
proximate property of the field equations (as usual for emergent systems), valid 
up to some order in perturbation theory. In particular, if we analyze the third 
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order contribution in (|129p we geiQ 



V {dXkf 



6 {dXkf 



• (140) 



This contribution is clearly not Lorcntz invariant if the solution ip pinpoints a 
specific direction, as for example when the Minkowski metric is emergent. As a 
matter of fact our theory will show aether like effects beyond second order. 

So far, we have hence generalized and extended results familiar to the ana- 
logue gravity community. However, as said, a typical drawback of analogue 
gravity models is related to the fact that they show only the emergence of 
a background Lorentzian geometry while they are unable to reproduce a ge- 
ometrodynamics of any sort. In what follows, we shall show that our model 
overcomes this drawback and indeed is able to describe the emergence of a the- 
ory for scalar gravity. This theory will come out to be the only known other 
theory of gravitation, apart from General Relativity, which satisfy the strong 
equivalence principle [69], i.e. Nordstrom gravity. 



6.2 Emergence of Nordstrom gravity 

In this section, we describe how we can recover a relativistic scalar gravity 
theory from a Lagrangian of the type (|126p . when ground state is such that 
the perturbations are living (at the lowest order in perturbation theory) in a 
Minkowski spacetime. So, let us start from the truncated Lagrangian for the 
perturbations (|139[) that we obtained in the previous section. This Lagrangian 
can simply be rewritten in terms of the (real) multiplet — (}pi, ^n) as 

C,s{^) = il^'''{d^^f{d,^). (141) 

This system has a global 0{N) symmetry which has emerged as well from the 
initial Lagrangian (jl26p . It is hence quite natural to rewrite the multiplet ip by 
introducing an amplitude characterized by a scalar field ^(x) and a multiplet 
4>{x) with N components such thal[f| 

f VI \ / 01 \ 

: =$ : , with |0|2 = ^02^£2^ (142) 

\ / \ 0Ar / 

i is an arbitrary length parameter to keep the dimension right. In particular, 
$ is dimensionless and (j) has the dimension of a length. $ is the field invariant 
under 0{N) transformations, whereas cj) does transform under 0{N). As we 
shall see, this field redefinition will provide us the means to identify gravity and 

^We are in the mono-metric case, so that F{Xi, X^) = f{Xi) + ... + f{X^), and 
ipk = i/", Vfc. 

^Our field redefinition is the generaUzation of the so-called Madelung representation 1121 . 



34 



matter degrees of freedom. The Lagrangian for the perturbations (|14ip reads 
now af0 

-Coff ("^1, ■■•Viv) — > -Ccff ($, 01, ...(/'Ar) = 

i 

where A is a Lagrange multipher. We recognize in particular the action for a 
non- linear sigma model given in terms of the fields (pi. The associated equations 
of motion are 

v'^'^ifd^d,'^ - $^a^0.a,0.) = 0, (144) 

i 

77^''(2a^$a,0, + <i>^d^d,^, + ^d^cPjd^P'^cj),) = 0, (145) 
|(/)p-£2 = 0. (146) 

If we introduce the (conformally flat) metric 

g^Ax) (147) 

the equations of motion p45p can be simply rewritten as 

Dgcl), + ■^ff^''a^0,a,0fc(5^V^ = 0, (148) 

where we have introduced the d'Alembertian Og for the metric g and used that 
\/—g = $^ and g^'^ = $"^77^''. Notice that equation (|145p can be rewritten in 
the form (|148p using the metric redefinition (|147p only in four dimensions. To 
be consistent, the change of variable $ g^i, should be completed with the 
constraint that g^j^ is conformally flat, that is 

Cc^pM = 0, (149) 

where Cap-yS is the Weyl tensor. 

Eq. (|148|) suggests that the gravitational degree of freedom should be encoded 
in the scalar field whereas matter should be encoded in the (pi. We are 
therefore aiming at a scalar theory of gravity with actions: 

Scff = j dx*y/^Ccff = S'grav + •S'matter, (150) 

^grav J rfa;^\/^r7'"5^*3,$, (151) 

^matter = J dx^ ^ ^^^'^ d^<P,d,cP, + X{\<P\^ - f)^ , (152) 



^We use the normalization condition |0p = I? , which implies in particular (j)id^(f>i = 0. 
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where we have exphcitly written the volume element ^J—t] = 1 so to make clear 
that these actions are given in flat spacetime. 

It is easy to see that the very same actions can be recast in the form of 
actions in a curved spacetime endowed with the metric (jl47p . In particular for 
the matter action in (|152p one has 



+\'m'-f)], (153) 



where we have suitably rescaled the Lagrange multiplier to A'. This allows to 
construct the stress-energy tensor r^,y for the non-linear sigma model, and its 
trace T with respect to the metric g: 

T^. = ^^^gJ^T = E {9^^^d.<l>, - ^.9,,.(g"'5aa0.5^0.)) , (154) 

i 

Finally, the above result, together with the recognition that the Ricci scalar R, 
associated to the metric g^v, can be written as R = — 6n,,$/$'^, allows us to 
rewrite Eq. (|144|) as the Einstein-Fokker equation 

i 

In summary, we can gather together the equations of motion (jl48p . (|149|) . 
(|155p . obtained by introducing the metric p47p . we have 

R =-^T, C^i3js = 0. (156) 

Dg(P, + ^g^''d^<f>jd,(f>kP''(t>, = 0, |0p -f = 0. (157) 
We recognize the equations of motion as those for Nordstrom gravity 

R = 247rGN T, Cap-yS = 0, (158) 

coupled to a non- linear sigma model. Indeed, the rewriting of (|144p - (|146p into 
the form (|156p - (|157p . is a special case of the procedure suggested by Einstein 
and Fokker so to cast Nordstrom gravity in a geometrical form [70] , 

We see from the above equation that the Newton constant Gn in our model 
has to be proportional to However, in identifying the exact relation between 
the two quantities, some care has to be given to the fact that the stress-energy 
tensors appearing respectively in equation (I156P and equation (llSSp do not share 
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the same dimensions. This is due to the fact that the fields have the dimension 
of a length rather than the usual one of an energy. This implies that in order 
to really compare the expressions one has to suitably rescale our fields with a 
dimensional factor, S, which in the end would combine with £ so to produce an 
energy, dim[€S] = energy. In particular, is easy to check that one has to assume 
47r£^S^ = -Bpianck i'^ order to recover the standard value of Gn (assuming c as 
the observed speed of signals and h as the quantum of action) . As a final remark, 
we should stress that the scale i is completely arbitrary within the emergent 
system and in principle should be derived from the physics of the "atoms of 
spacetime" whose large N limit gives rise to (|126p . 

Accidentally, the above discussion also shows that, once the fields are suit- 
ably rescaled so to have the right dimensions, the constraint appearing in 
Ea. (|157p is fixing the norm of the multiplet to be equal to the square of the 
Planck energy. This implies that the interaction terms in the aforementioned 
equation are indeed Planck-suppressed and hence negligible at low energy. This 
should not be a surprise, given that in the end £S is the only energy scale present 
in our model. It is conceivable that more complicate frameworks, possibly en- 
dowed with many dimensional constants, will introduce a hierarchy of energy 
scales and hence break the degeneracy between the scale of gravity and the scale 
of matter interactions. 

6.3 Lessons from the toy model 

Initially, we have considered fields that live in a Euclidean space, and showed 
that there exists a class of Lagrangians such that the perturbations around some 
classical solutions "0 propagate in a Minkowski spacetime. In this case ip is essen- 
tially picking up a preferred direction, so that we have a spontaneous symmetry 
breaking of the Euclidean symmetry. The apparent change of signature is free 
of the problems usually met in signature change frameworks since the theory is 
fundamentally Euclidean. Lorcntz symmetry is only approximate, and in this 
sense it is emergent. 

The main lesson we want to emphasize here is that Lorentzian signature can 
emerge from a fundamental Euclidean theory and this process can in principle be 
reconstructed by observers living in the emergent system. In fact, while from 
the perturbations point of view it is a priori difficult to see the fundamental 
Euclidean nature of the world, this could be guessed from the fact that some 
Lorentz symmetry breaking would appear at high energy (in our case in the 
form of a non-dynamical ether field). 

In the second part of the discussion of this toy model, using a natural field 
redefinition adapted to the symmetries of the system, we have identified from the 
perturbations ipi , a scalar field $ encoding gravitational degrees of freedom and 
a set of scalar fields 0i (a non-linear sigma model) encoding matter fields. In this 
sense, gravity and matter are both emergent at the same level. This approach 
is then rather different from the one of analogue models of gravity where one 
usually identifies the analogue of the gravitational degrees of freedom with the 
"background" fields, i.e. the condensate or the solution -0 of the equations 
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of motion. Indeed, following this line of thought in looking for a theory of 
gravitational dynamics, we would be led to require that the fundamental field 
theory (|126p must be endowed with diffeomorphisms invariance from the very 
start — the symmetries of the background are identical by construction to the 
ones of the fundamental theory. This would imply that one would have to obtain 
gravity from a theory which is already diffeomorphisms invariant and hence 
most probably with a form very close to some known theory of gravitation. For 
these reasons, we do expect that if an emergent picture is indeed appropriate 
for gravitation, then it should be of the sort presented here, with both matter 
and gravity emerging at the same level. Of course, it is not possible to exclude 
that a full fledged theory of gravity could emerge, together with the notion of 
manifold, in a single step from the eventual semiclassical/large number limit 
of the fundamental objects. In this case, however, we would still have a very 
different picture from the one envisaged in analogue models of gravity. 

In particular, this allows not only for an emergent local Lorentz invariance for 
the perturbations dynamics but it leads as well to an emergent diffeomorphisms 
invariance. In fact, we saw how the equations of motion (|144p and (|145p could be 
rewritten in a completely equivalent way using a conformally flat metric (|147[) . 
Most noticeably, they can be rewritten in an evidently diffeomorphisms invariant 
form, from the point of view of "matter fields observers" . In fact, following the 
standard hole argument (see [7T] for a careful discussion of the various issues 
related to diffcomorphism invariance), this also implies that the coordinates 
x'^, used to parameterized our theory, do not have any physical meaning from 
the point of view of the (pi "matter observers" . They are merely parameters. 
In agreement with the fact that diffeomorphisms invariance is emergent in our 
system, it can be noted that the cubic contribution ()140p ends up breaking it 
at the same level it breaks Lorentz invariance. 

Furthermore, Nordstrom gravity is also a nice framework for discussing the 
subtle distinction between background independence and diffeomorphisms in- 
variance [75] . We call background some geometrical degrees of freedom that are 
not dynamical. For example, in General Relativity the topology of the manifold 
and its dimension, or the signature of the metric, can be considered as (triv- 
ial) background quantities. We can therefore have some specific background 
structures while still having diffeomorphisms invariance. Nordstrom gravity is 
encoded in conformally flat metrics. If one considers fields which are confor- 
mally coupled to the metric (such as the electromagnetic field), these fields only 
see the metric 77^;^ which is of course not dynamical. The Minkowski metric 
can be see then as a background structure, this is what one may call a "prior 
geometry" (e.g. see [30]). One may hence say that diffcomorphism invariance is 
somewhat of a weaker form in Nordstrom gravity with respect the one present 
in general relativity. 

In particular, while the essence of diffcomorphism invariance in General Rel- 
ativity is encoded in the associated Hamiltonian constraints, these are not de- 
fined in the present formulation of Nordstrom gravity. Furthermore, in the most 
general implementations of Norstrom theory, quantities can be built which man- 
ifestly include the background structure rj^i, and hence are not diffcomorphism 
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invariant. However, within our model, the prior geometry cannot be detected. 
Indeed, in order to detect the Minkowski background, one should be able to pro- 
pose a method to pinpoint the conformal factor $^ in the relation g^^ = ^^r^^y. 
However, a careful analysis shows that this is actually impossible. Let us elab- 
orate on this point. If we perform a conformal transformation, — > x^{x), the 
equations of motions associated to (|14ip are transforming like 

□^V^, = -> Uf^if^ = 0, (159) 

where rj and ry are two different Minkowski metrics related by some conformal 
factor \{x). Therefore, rj and fj are indistinguishable, due to conformal invari- 
ance the equations of motion for ipi. Hence, what appears to be a background 
structure, namely "q^vi is ambiguously defined, and the coordinates in which 
the equations of motion for the fields Lpi are written have no operational mean- 
ing, they are mere labels. Furthermore, this ambiguity in the definition of what 
would be called a background structure implies an ambiguity on the definition 
of the conformal factor relating the physical metric to the would-be background 
structure. In this sense, within this very specific implementation of the model 
which has conformal invariance, there is no Minkowski geometry as a back- 
ground. There is a background structure, which is the conformal structure of 
Minkowski spacetime. This is a mild limitation of our simple toy model as a 
diffeomorphism invariant, background independent system. 

Of course, the above discussion holds only at the lowest order in the fields 
Lpi. As previously discussed, higher orders in perturbation theory will generate 
terms like (jl40p producing a breaking of the conformal symmetry and hence the 
appearance of the background structures, i.e. the Euclidean space and the d^iil) 
which have selected the timelike direction. 

Finally, Nordstrom gravity is only a scalar gravity theory, which has been 
falsified by experiments (e.g. the theory does not predict the bending of light). 
In order to obtain a more physical theory, in particular General Relativity, 
one should surely look for more complicated emergent Lagrangians than (|126[) . 
Of course, one would in this case aim to obtain the emergence of a theory 
characterized by spin-2 gravitons (while in Nordstrom theory the graviton is 
just a scalar). This would open a door to a possible conflict with the Weinberg- 
Witten theorem (see section [01 and [36 J- However, there are many ways in 
which such a theorem can be evaded and in particular one may guess that 
analogue models inspired mechanisms like the one discussed here will generically 
lead to Lagrangian which show Lorentz and diffeomorphism invariance only as 
approximate symmetries for the lowest order in the perturbative expansion. 

It is unclear which sort of generalization may still lead to some viable grav- 
itational theory from the perturbations dynamics. For example, the simple 
addition of a potential will in general prevent the selection of a preferred direc- 
tion, except in regions where the potential is almost flat. Moreover, it would 
also spoil the metric interpretation of the theory. For example, the terms |(^|" 
for n > 1 and ^ 4 cannot be rewritten as an interaction between the matter 
field fields living on the conformal metric ^^rj^^, when using the change of 
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variables (|142p (although it is interesting to note that a \ip\'^ term would give 
Nordstrom gravity with a cosmological constant). 

Perhaps, this toy model is too simple: in order to obtain more realistic 
theories one has to use more ingredients. Nevertheless, the reader should be 
convinced, by now, that such an objective could be not completely hopeless, as 
this primitive toy model suggests. Of course, much more should be done in this 
direction. 

6.4 The role of symmetries 

It is interesting to discuss in details the features that allowed the construction of 
such a toy model. In particular, it is important to stress the role of symmetries, 
in order to make clear the way in which they enter at the various levels. As 
in the case of selecting Riemannian geometry out of Finsler geometry, here 
there are some symmetries which are absolutely essential: it is only due to 
their presence that we do have an emergent gravitational system possessing a 
geometrical nature. 

We have seen that in order to produce a working model, a number of prop- 
erties must be assumed. First of all, there is an underlying /S'0(4) symme- 
try which allows us to use particularly simple afhne solutions. This IS0{4), 
when spontaneously broken, can lead to an approximate Poincare invariance. 
Moreover, the masslessness of the resulting modes is promoting this Poincare 
invariance to a full conformal invariance, which is approximate as well. This 
conformal invariance is the key symmetry which hides the background struc- 
ture, forbidding a low energy observer to detect a background metric structure 
(there is only a background conformal structure). 

Conformal invariance seems to be deeply intertwined with the possibility 
of writing down the resulting equations of motion in the form of a system of 
diffeomorphism invariant equations, as we have seen. However, in order for the 
Lagrangian (|139p to be conformal invariant, there must be an overall 0{N) 
symmetry between the fields. This symmetry is just the other side of the coin 
of the mechanism leading to the monometricity. If two fields move in different 
metrics, clearly this 0{N) is broken and the entire model fails to provide a 
geometric picture, let alone a diffeo-invariant one. 

In general, one should expect that in any situation in which the metric is 
an emergent structure, there should be a mechanism taking care of the fact 
that different matter fields propagate over the same geometry. In this picture, 
where a manifold is given from the beginning, the role of internal and spacetime 
symmetries is crucial. The behavior we have described is not general at all. Of 
course, one could conclude that this kind of models is somehow contrived and 
unnatural. 

However, there is also a positive side: given that symmetries (both of the 
equations of motion and of the ground state) play a crucial role in the emergence 
mechanism, the fact that our universe seems to be ruled, at large scales, by 
general relativity and locally by special relativity, suggests that not all the pre- 
geometric scenarios are viable, and that there are rather strong constraints on 
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what are the possible mechanism of emergence. 



7 Conclusions 

In summary, we hope that the two examples of emergent gravitational dynamics 
presented here have suitably illustrated the potentialities of emergent gravity 
models inspired by the analogue gravity perspective. The first case, the one of 
a BEC system, has shown us that a gravitational-like dynamics (with a small 
cosmological constant) seems to be a natural by product of a condensation mech- 
anism. However, the analogy with the real world was not only limited by the 
Lorentz breaking scale but also from the fact that gravity and matter seem to be 
living at rather different levels (gravity is the condensate while matter is asso- 
ciated to the quasi-particle states). The toy model discussed in the second part 
of this proceedings is aimed at overcoming this problems. It does have Lorentz 
violation but only in the limit of large fluctuations of the fundamental fields (jji- 
This implies that the gravitational dynamics is no more endowed with a massive 
graviton whose mass scale is set by the UV Lorentz breaking scale of the system. 
Furthermore, this toy models shows how time and diffeomorphism invariance 
might emerge. In particular the latter is allowed by the special symmetries of 
the system. As we have said previously, similar symmetries are probably needed 
anyway to protect the IR limit of the theory from large violations of Lorentz 
invariance of the equations. Given that the presence of such violations seems 
to be a very natural way around to the Weinberg- Witten theorem obstruction, 
it might be that the next step towards a satisfactory emergent gravity scenario 
might have to consist in finding which sort of mechanism, possibly a custodial 
symmetry, could simultaneously guarantee the background independence of the 
emerging dynamics as well as a very accurate Local Lorentz invariance of the 
emergent spacetime. We hope to address these questions in future work. 
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